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ABSTRACT 

We analyze high resolution, N-body hydrodynamical simulations of fiducial galaxy clusters to probe 
tidal stripping of the dark matter subhalos. These simulations include a prescription for star formation 
allowing us to track the fate of the stellar component as well. We investigate the effect of tidal 
stripping on cluster galaxies hosted in these dark matter subhalos as a function of projected cluster- 
centric radius. To quantify the extent of the dark matter halos of cluster galaxies, we introduce the 
half mass radius ri/2as a diagnostic, and study its evolution with projected cluster-centric distance 
Ras a function of redshift. We find a well defined trend for {ri/2,R) '■ the closer the galaxies arc to 
the center of the cluster, the smaller the half mass radius. Interestingly, this trend is inferred in all 
redshift frames examined in this work ranging from z = to z = 0.7. Atz = 0, galaxy halos in the 
central regions of clusters are found to be highly truncated, with the most compact half mass radius 
of lOkpc. We also find that ri/2 depends on luminosity and we present scaling relations of ri/2 with 
galaxy luminosity. The corresponding total mass of the cluster galaxies is also found to increase 
with projected cluster-centric distance and luminosity, but with more scatter than the (?'i/2, R) trend. 
Comparing the distribution of stellar mass to total mass for cluster galaxies, we find that the dark 
matter component is preferentially stripped, whereas the stellar component is much less affected by 
tidal forces. We compare these results with galaxy-galaxy lensing probes of and find qualitative 
agreement. Future surveys with space based telescopes such as DUNE and SNAP, that combine wide 
field and high resolution imaging, will be able to probe the predicted (ri/2, R) relation observationally. 

Subject headings: galaxies: dark matter halos - numerical simulations: N-body, hydrodynamical 



1. INTRODUCTION 



The dependence of galaxy properti es on 
ment is well established (see 



en viron- 
Adami et all 119981: 



Lanzoni et all 120051 : iBoselh fc Gavazzil 120061: iDe Lucial 



20061) . One of the most extreme environments for galax- 



ies is inside a massive galaxy cluster, where active, strong 
tidal forces are exerted by the global cluster potential. 
The theoretical expectation is that the global tidal field 
of a massive, dense cluster potential is strong enough 
to truncate the dark matter halos of gal axies that tra- 
verse the cluster core. Early work (e.g. [Merrittl Il983t 
iMalumuth fc Richstong 119841 : iMerrittI Il984f ) found that 
a large fraction of the mass initially attached to galax- 
ies in t he central Mpc is s tr ipped. lAvila- Reese et al.l 
dHH); iBullock et all ((200ll) : lAvila-Reese et al.l (|2005[ ) 



found that halos in dense environments are more trun- 
cated and more compact than their isolated counter- 
parts of the same luminosity. On galaxy scal e s, us - 
ing dark matter only simulations. iDiemand et al] (|2007l ) 
have studied the evolution of subhalos in the Via Lactea 
host halo, and find that tidal forces remove subhalo mass 
from the outside in, which leads to higher concentrations 
for subhalos located in the inner regions compared to 
field halos of the same mass. Detailed studies of the evo- 
lution of subhalos in clusters i n dark matter only s imula- 
tions have been performed bv lGhigna et al.l (jl998L here- 
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after G98) and iDe Lucia etal\ (j2004l ). G98 find that 
the dominant interactions in cluster environments are be- 
tween the global cluster tidal field and individual galaxies 
after z — 2, and that the cluster tidal field significantly 
strips galaxy halos. Moreover, both nu merical simula- 
tions (G98) and analytical calculations (|Mamonl 120001 ) 
predict that the tidal radius of a given galaxy depends 
on its cluster-centric distance. As a consequence, the 
closer to the center the galaxies are, the stronger are 
the tidal forces they will experience, resulting in more 
compact subhalos. G98 probed the characteristic extent 
of galaxy subhalos with cluster-centric distance, at z = 
and z = Q.5. Considering the three dimensional cluster- 
centric distance, they find that subhalo extents decreases 
towards the cluster center, but this trend was weak and 
therefore hard to detect at z = 0.5. They argue that this 
is due to the fact that at 2: = 0.5, the cluster has quite 
an anisotropic mass distribution and tides are efficient 
only at its very center. When considering the projected 
cluster-centric distance, this trend does not totally dis- 
appear but becomes marginal, suggesting that subhalo 
extents do not strongly depend on present day projected 
cluster-centric distances. 

From an observational point of view, gravitational 
lensing seems the only currently viable method to 
probe the extent of dark matter subhalos in clusters, 
and this statement is specially true for cluster galax- 
ies. In the case of field gal axies, satell i te dy nam- 
ics can be used as well {e.g. iPrada et al.l l2003f ) and 
gives results consistent with current galaxy-galaxy lens- 
ing studies. The deflection caused by galaxy scale 
mass concentrations is small (quantified by a shear 
7 ~ O.Of) and thus is challenging to detect. However, 
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galaxy- galaxy lensing studies in clusters h ave bee n per- 
formed iNataraian et a l.l [1998l:lGeiger fc S clincidcj|1999t 
INataraian et al. l2002al lb : iLimousin et al.l i2007a) where 



this shear signal is boosted by the smoothly distributed 
large scale dark matter distribution. These studies have 
successfully statistically detected the weak lensing signal 
generated by clust er galaxies (see also recent work by 
iHalkola et al.l [20071 ). These analyses have provided evi- 
dence for truncation of galaxy dark matter halos in high 
density environments. The inferred typical half mass ra- 
dius was found to be typically more compact than 50 kpc, 
whereas half mass radii larger than 200 k pc are derived 
for fi e ld galaxies of equiva l ent lum inosity (iFischer et al.l 
[2OOOI: iMcKav et all [2OOII: [Ho^stra et aP 120031 I2004D . 
Moreover, when using the NFW (Na varro et al.lll997f ) pro- 
file, co ncentration parameter s greater than 20 were in- 
ferred (jLimousin "etall [2OO6I) for some clusters. These 
results are in qualitative agreement with the tidal strip- 
ing scenario. Recently, using a large hst mosaic cover- 
ing up to 5 Mpc from the center of galaxy cluster C10024 
at z = 0.39, Natarajan et al. (2008) were able to probe 
the galaxy population in three radial bins and inferred a 
larger extent for the halos of galaxies living in the out- 
skirts of the cluster (i.e. at a cluster-centric distance 
between 3 and 5 Mpc) compared to the galaxies living in 
the core of the cluster (here between and 3 Mpc) . 

In this paper, we want to investigate how tidal forces 
shape the fate of galaxy dark matter halos, and in par- 
ticular, we focus on the extent of cluster galaxy dark 
matter halos and probe its evolution with cluster-centric 
distance and redshift. We analyze N-body hydrodynam- 
ical simulations of two fiducial galaxy clusters that con- 
tain not only dark matter particles but also stars and gas 
particles that interact through many physical processes. 
In order to make predictions for future surveys that will 
constrain the extent of the halos as a function of cluster- 
centric distance, we consider the projected cluster-centric 
distances since this is what can be probed observation- 
ally. The extent of subhalos is quantified by the half mass 
radius ri/2 since it is directly comparable to the charac- 
teristic extent probed via galaxy-galaxy lensing studies. 

This paper is organized as follows: in Section 2, we 
discuss the properties of the two simulated clusters stud- 
ied here. In Section 3, we investigate how the half mass 
radius r]^/2Etnd the total mass Mtot evolve with cluster- 
centric distance, for different assembly stages (redshift 
frames) in these simulations. In Section 4, we investi- 
gate how ri/2and Mtot evolve with luminosity. We fit 
scaling relations for ri/2 , Mtot and the luminosity, and 
present total mass to light ratios. In Section 5, we in- 
vestigate the ratio of the total to the stellar mass and 
its evolution with cluster-centric distance and redshift. 
In Section 6, we compare galaxy-galaxy lensing results 
to our theoretical predictions, and finally we present our 
conclusions in Section 7. All our results in this paper are 
scaled to the flat, low matter density Acdm cosmology 
with = 0.3, r^A — 0.7 and a Hubble constant Ho = 70 
kms~^ Mpc~^. 

2. SIMULATIONS OF OUR FIDUCIAL CLUSTERS 

2.1. Cluster Properties 

N-body hydrodynamical (treesph) simulations of the 
formation and evolution of two galaxy clusters in a Acdm 



cosmology have been performed. The simulations in- 
clude metallicity-dependent radiative cooling, star for- 
mation according to different initial mass functions, en- 
ergy feedback as strong starburst-driven galactic super- 
winds, chemical evolution with non-instantaneous recy- 
cling of gas and heavy elements, effects of a metagalactic 
ultraviolet field and thermal conduction in the intraclus- 
ter medium. For full details on these simulations and for 
a comparison of the properties of the simulated clusters 
to obse rvations, we refer the rea d er to a series of three 
paper s: iRomeo et all (l200l[2005h :[S ommer-Larsen et al.l 
(|2005( 1. hereafter RSL. 

The main difference between the two simulated clus- 
ters is their temperature (or equivalently their mass) and 
the numerical resolution. The more massive cluster is a 
6keV cluster, that we will hereafter refer to as COMA. 
At z — 0, its virial radius is equal to 2.9 Mpc, and the 
corresponding virial mass is equal to 1.3 x 10^^ Mq. The 
second cluster is a 3.1 keV cluster that we will refer to as 
VIRGO. At z =0 , its virial radius is equal to 1.8 Mpc, and 
the corresponding virial mass is equal to 2.8 x 10^^ M©. 
The particle mass in viRGO is 8 times lower than the one 
used to simulate COMA . Thus the VIRGO simulation is a 
higher resolution one and we are able to resolve smaller 
galaxies. On the other hand, since COMA is a more mas- 
sive cluster, more massive galaxies than in viRGO will 
form. 

The highest redshift frame we analyze is z = 0.7. At 
this redshift, both clusters are virialized and there is a 
clear dominant cluster halo. At higher redshifts, the clus- 
ters are not well defined (no clear dominant halo) and 
are still under construction with major merger events. 
However, for any redshift slice, COMA is less relaxed than 
VIRGO. In the case of COMA at z = 0.7, simulations re- 
veal some infalling galaxies, whereas that is not the case 
in VIRGO . This difference in dynamical states can be 
quantified by the redshift of the last major merger event: 
for COMA , the last merger event occurs at around z = 0.8, 
whereas it is at around z = 1.5 for viRGO . This means 
that COMA is dynamically younger than viRGO , and, as 
we will see, the viRGO galaxy population is more relaxed, 
with less scatter in the relations we study in this work. 
Due to the difference in dynamical states between the 
two simulated clusters, it is not easy to directly compare 
their respective properties, in the sense that COMA is not 
simply a rescaling of viRGO at a higher temperature. Pic- 
tures and movies of the simulated clusters can be found 
at http://www.tac. dk~jslarsen/Clusters/ 

2.2. Extracting subhalos from the simulated clusters 

The first step of our analysis is to extract substruc- 
tures from the simulated cluster in order to get a catalog 
of cluster galaxies. This procedure is detailed in RSL and 
summarized below. In order to ascertain that we iden- 
tify all galaxies, and hence sub-halos, in the main cluster 
halo the following approach was adopted. A cubic grid 
of cube-length A^=10 kpc is overlaid the cluster, and all 
cubes containing at least A^th=2 star particles are identi- 
fied. Subsequently, each selected cube is embedded in a 
larger cube of cube-length 3AZ. If this larger cube con- 
tains at least 7Vinin=7 star particles, which are gravita- 
tionally bound by its content of gas, stars and dark mat- 
ter the system is identified as a potential galaxy. Since 
the method can return several, almost identical versions 
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of the same galaxy only the one containing the largest 
number of star particles is kept and classified as a galaxy. 
We tested the galaxy identification algorithm by varying 
Al, Nth-, -^min, and also the numerical resolution, and 
found it to be adequately robust for the purposes of this 
paper. 

2.3. Suhhalos Properties 

Once the subhalos have been identified, galaxy half 
mass radii, ri/2 , and total masses, Mtot , were deter- 
mined by an iterative approach, described below. For 
a given galaxy the density of the smooth cluster "back- 
ground" of dark matter and hot, intra-cluster gas at the 
position of the galaxy is determined. This is done by cal- 
culating the average total density in 5 kpc thick shells of 
radii from to 1.5 Tj. The "search radius" , r-g, is initially 
set to Art, where rj, is the tidal ra dius of the galaxy at 
its position in the cluster defined by (|Binnev fc Tremaind 
[1988,): 

All shells containing other substructures, i.e., galax- 
ies and/or dark matter halos, are subsequently removed 
using a background density fluctuation criterion, and 
finally the background density is determined from the 
remaining shells. Second, the cumulative mass of the 
galaxy, including its dark matter halo (above the back- 
ground density) is determined in 2.5 kpc thick shells, 
going from inside and out, until a shell is reached of den- 
sity equal to or less than that of the background. Using 
the cumulative mass distribution, the half mass radius 
ri/2 , is subsequently determined by linear interpolation. 
Third, the search radius is set to Ari/2 , and the above 
procedure is redone. This is done repeatedly until con- 
vergence. Tests with different parameter values show the 
above procedure to be an efficient and robust way for 
determining half mass radii and total masses of cluster 
galaxies. 

By construction, r 1/2 is the radius within which the 
3D mass of the simulated galaxies equals half of the to- 
tal mass. We will refer to this quantity as being the 
3D half mass radius. We also pursue a 2D analysis. In 
particular, we looked at R1/2 , the radius within which 
the 2D mass equals half of the total mass, as described 
below. Once we know the 3D region (assumed spheri- 
cal) where a given galaxy is confined and contained, we 
projected the excess mass density in this region (above 
the cluster mean density) onto the plane along the three 
cardinal directions, and averaged over the results. This 
gives us a robust estimate of the 2D projected mass, that 
we use to calculate R1/2 ■ We will refer to this quantity 
as the 2D half mass radius. We will compare these quan- 
tities in Section 3. However, when looking at possible 
scalings of the half mass radius with luminosity, we will 
use rx/2 0iily. The quantity ri/2 ''^^ robust as it is easily 
calculated and well defined for most popularly used mass 
models. Moreover, ri/2 is not the result of an average of 
projected quantities, thus in principle a less noisy quan- 
tity than Ri/2 . 

There is a large scatter in the properties of the galaxies 
that form in these cluster simulations. This scatter is 
mainly due to the diverse orbital histories of galaxies. To 
consolidate results, we define bins, in which we compute 



the median of the points as well as the standard deviation 
a. The error is defined as a/\/^, where N is the number 
of objects in a given bin. In the following, when needed, 
we fit a linear function to the points corresponding to 
the galaxies, and not to the binned data. 

2.4. Projected cluster- centric distance 

Since observational probes are sensitive to projected 
cluster-centric distances, we project the 3D radial dis- 
tance of the galaxies with respect to the center of the 
cluster over the three cardinal axes, and our cluster 
galaxy catalog is made by merging the three catalogs 
corresponding to each projection. We thus artificially 
increase the data by a factor of three. This is not unrea- 
sonable to do since we can imagine that three different 
observers have been looking at the cluster from the three 
cardinal axes and have merged their corresponding data 
sets. This procedure is equivalent to observing three dif- 
ferent clusters in a similar dynamical state. It is impor- 
tant at this point to make sure that the relation between 
half mass radius and projected cluster-centric distance 
we aim to probe does not depend on the adopted view- 
ing angle and that we do not introduce any systematics 
by merging catalogs corresponding to each projection. 
To this purpose, we have looked at the evolution of the 
half mass radius with respect to the three different pro- 
jections of the cluster-centric distance. We find no bias. 
To illustrate this point, we show on Fig. [l]the evolution 
of the half mass radius as a function of the projected 
cluster-centric distances, considering the binned data and 
the VIRGO cluster. We see that the relation does not de- 
pend on the adopted projected cluster-centric distance. 
The same conclusion can be drawn when considering the 
unbinned data points. Moreover, we reach similar con- 
clusions in the case of the COMA cluster. Also shown for 
comparison is the relation when merging the data sets 
as done in the rest of the paper. Note that looking at 
possible variations with viewing angle is interesting from 
an observational perspective and suggests that projection 
effects are not important for observational studies. 

3. TRENDS WITH CLUSTER-CENTRIC DISTANCE 

The extent and efficiency of tidal stripping depends pri- 
marily on the details of the orbital parameters of the sub- 
halo. In this Section, we study the evolution of the half 
mass radii {ri/2 , R-1/2 ) and the total mass Mtot as a func- 
tion of the projected cluster-centric distance R, for dif- 
ferent redshift frames. In Fig. [21 we plot ri/2 and R1/2 as 
a function of R both for viRGO and COMA , for the two 
redshift frames we consider, i.e. z = and z = 0.7. We 
clearly see that both the (ri/2,R)and (R1/2, R) trends 
are well defined, at any redshift. Closer to the center of 
the cluster the extent of the dark matter halos of the clus- 
ter galaxies is smaller. Thus closer to the center of the 
cluster, stronger tidal stripping is experienced by galax- 
ies. The (ri/2jR)and (R1/21 R) trends are illustrated on 
Fig. [2] for the two extreme redshift frames considered in 
the simulations, z = Q and z = 0.7. This choice is moti- 
vated by clarity of the figure and to enhance the evolution 
between z = and z = 0.7. These trends are inferred in 
all redshifts frames extracted from the simulations (z ~ 0, 
z = 0.2, z = 0.4 and z — 0.7). Note that if this trend is in- 
ferred both at z = and z — 0.7, then it must be present 
for intermediate redshifts as well. 
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Fig. 1. — median half-mass radii r^j-^ as a function of projected 
cluster-centric distance in VIRGO . This cluster-centric distance cor- 
responds to the projection of the 3D cluster-centric distance along 
the three cardinal axes: along the x axis (green), the y axis (red) 
and the z axis (blue). Also show is the relation when merging the 
three data sets (black) which corresponds to the solid line in the 
upper left panel of Fig. 2. We see no bias with viewing angle. 

We see that for any given R , r-^j^ and R1/2 are system- 
atically smaller at lower redshift compared to the high 
redshift frame: by z = 0, the galaxies have had more pas- 
sages through the cluster center, have been stripped more 
and thus have a smaller extent. As expected, we see that 
the 2D half mass radius Ri/2is systematically smaller 
than the 3D half mass radius r^j^ ■ 

In Fig. [21 we plot the total mass Mtot as a function of 
R both for viRGO and COMA , for the two redshift frames 
we consider. At z = 0, we infer a (Mtot, R) trend, but 
the scatter is larger than the one associated with the 
(ri/2, R) trend. This difference in trends can be under- 
stood as follows: consider a galaxy of mass m with a 
circular velocity Vc located in the cluster whose mass is 
M, we can relate the scatter on the tidal radius {6rt/rt) 
to the scatter on the galaxy mass {5m /m). The tidal 
radius of the galaxy is proportional to: 

(Eq. [T|) and the mass of the galaxy is given by: 

m (X rt vl, (3) 

thus Tt oc vV^ , so rt oc Vc and 5rt oc 5vc- 
Differentiating Eq. [31 we can write: 

5m ^ 5rt + rt 2 Vc 5vc oc 3 5vc vl (4) 

Combining Eq. [3]and[31 we get: 

5m Z5vcvl 5rt 

— = — = 3 — 5) 

m rt rt 

Therefore finding a larger relative scatter in the (Mtot , 
R) trend than in the (ri/2, R) trend is to be expected. 
Note that at z = 0.7, coma is not well relaxed and some 
galaxies have not experienced the tidal field of the clus- 
ter potential yet. As a consequence, the outer radial 
bin contains some infalling galaxies which have values of 



'"1/2 and total masses typical of field galaxy values (i.e. 
ri/2 > 200 kpc and Mtot > 10^^ Mq as obtained from 
galaxy-galaxy lensing and satellite dynamics analysis, see 
Section 1). 

In Fig. 2, we also show the results from G98 at z = 0. 
G98 studied a simulated galaxy cluster with virial radius 
1.95 Mpc and virial mass 2.8 x 10^'' M0. These proper- 
ties makes it very similar to our viRGO cluster. Thus we 
compare their results (data points with triangle) to ours 
on VIRGO , for the redshift frame at 2: = (solid line). We 
find significant differences between our results and the 
one from G98. Before discussing these differences, we 
want to caution that both works are different and diffi- 
cult to compare. Therefore, the disagreement between 
these studies cannot be used to assess quantitatively the 
impact of the baryonic component in the tidal stripping 
process. We also note that there could be significant dif- 
ferences in subhalo finding algorithms. Kepping that in 
mind, three statements can be made from this compar- 
ison: data points from G98 have larger error bars and 
are systematically higher than ours. Moreover, the clear 
trend we infer appears marginal in the work by G98. 
The difference in the error bars comes from the statis- 
tics: within 2 Mpc, G98 identified 200 halos whereas 
we have ^ 650 halos within the same radius. The dif- 
ference in the value of ri /2 can be partly due to the dif- 
ference in resolution: G98 use particles whose mass is 
8.6 X 10^ Mq whereas we use particles whose masses are 
4.4 X 10^ Mq and 3.2 x 10^ M© for the stellar and dark 
matter respectively. The mass of our individual parti- 
cles is 3-20 times smaller, thus we resolve smaller galax- 
ies than G98: they adopt a minimum of 32 particles to 
define a halo whose individual properties are relevant, 
which translates into a minimum mass of 2.7 x lO"'^" M©. 
In VIRGO we have 477 galaxies of total mass smaller than 
this threshold. Therefore our sample is more complete 
in the inclusion of lower galaxy masses compared to the 
G98 sample. Since smaller galaxies tend to have smaller 
extents, this could explain the shift between our results. 

In order to properly investigate the role of the bary- 
onic component in the tidal stripping process, we would 
need a devoted dark matter only simulation of viRGO and 
COMA . Such simulations are underway, and will be pre- 
sented and analysed in a forthcoming publication, where 
we will be able to quantify the differences expected. 

4. SCALING OF ri/2 AND Mtot WITH LUMINOSITY & 
TOTAL MASS TO LIGHT RATIO 

Para metric strong l ensing s tudies (e.g. iKneib et al.l 
19961: lTvson etal.1 119981 : iNataraian et all Il998t 
Limousin et al I I2007bn do include the galaxies liv- 
ing in the core of the cluster in the modeling since they 
can locally pertub strong lensing features. In order to 
reduce the number of free parameters in lens modeling, 
the standard approach consists of scaling the galaxy 
parameters with luminosity as follows: 

ri/2 cx & Mtot oc (6) 

Galaxy-galaxy lensing studies also use similar scaling 
laws. The scaling usually used for the half mass radius 
in lensing studies is: a = 0.5, which assumes that the 
mass to light ratio is constant for all galaxies. At the 
present time strong lensing studies are not in a posi- 
tion to discriminate between different scaling relations 
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Fig. 2. — Upper panels: median half-mass radius (3D) and R1/2 (2D) as a function of projected cluster-centric distance, for 2 = and 
2 = 0.7; VIRGO (left) and COMA(right). For both clusters, the (riy2! R-) S'lid (R1/21 R) trends arc clearly defined in both redshift slices. Also 
shown for comparison with VIRGO are the data point from G98 at 2 = 0. As wc could have expected, we see that the 2D half mass radius 
1^1/2 is systematically smaller than the 3D half mass radius . Lower panels: median total mass as a function of projected cluster-centric 
distance, for 2 = and z = 0.7; virgo (left) and coma (right). For both clusters at 2; = 0, we find a well defined (Mtot , R) trend, though the 
scatter is larger than in the (ri/2>R)or (R1/2 > R) trends. At high redshift, the scatter is larger and the (Mtot, R) trend becomes marginal. 
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(|Halkola et al.ll2007t [Jullo et al.ll2007[ ). In order to study 
the possible value of the exponents defined by Eqn. [6] 
that could be used in strong lensing modeling, we want 
to examine the evolution of the half mass radius ri /2 and 
the total mass Mtot with luminosity L, as a function of 
redshift. We consider the galaxies located in the core 
of the cluster, (i.e. the galaxies that satisfy R< 500 
kpc) since strong lensing modeling deals with the inner 
part of a galaxy cluster (typically the central 100"). In 
Fig. [3] (upper panels), we plot ri/2 as a function of lumi- 
nosity L both for viRGO and COMA , at z = and z — 0.7. 
We fit a linear relation to the unbinned data points for 
each redshift frame (i.e. z = 0, 0.2, 0.4, 0.7) in order to 
probe the exponent of the scaling relations defined by 
Eqn. [6l The results are presented in Table [1] and the 
best linear fit is plotted on Fig. [31 We also fit a linear re- 
lation to the unbinned data points corresponding to the 
(Mtot, L) trend (Table [TJ. We find increasing (ri/2 , L) 
and (Mtot, L) trends at both redshifts, for the COMA and 
VIRGO central populations. Looking at the evolution of 
a and 5 with redshift, we do not infer any clear trend of 
evolution of these values with redshift, thus we average 
the exponents from the different redshift frames and we 
consider the averaged values in what follows. Comparing 
VIRGO and COMA , we see that the scatter is smaller in 
the case of virgo which is more relaxed than COMA (see 
error bars on the linear fitting results in Table [T|) . 

We caution that we do not claim to have derived fun- 
damental scaling relations here, but rather have fitted to 
our data sets. As can be seen in Fig. [3l it is clear that 
the relations likely have a more complicated form than a 
linear fit. 

These inferred scalings can be compared to the fun- 
damental plane derived for early-type galaxies in clus- 
ters. The values of 5 derived in this work are in good 
agreement with the fundamental plane analysis at z = 
for ea rly-type galaxies {5 — 1.35 ± 0.15: iJorgensen et al.l 
[l99fih . 

On Fig. [3] we plot the total mass to light ratio 
Mtot/L as a function of luminosity. We find that 
Mtot/L can reach values as high as 55, with a mean value 
around ^ 20-30. It is interesting to compare these total 
mass to light ratios to stellar mass to li ght ratios (M/L) 
inferred ob servationally. Recent studies (jJorgensen et al.l 
l2006l l2007l ) present a relation between the galaxy masses 
and the mass to light ratios, for galaxy clusters in the lo- 
cal Universe {z = 0) and in the high r edshift Universe 
(z = .8 - 0.9). The relations given bv iJorgensen et aP 
were derived in the B band. We converted their 
B band luminosities into the R band luminosities using 
a constant (b - R = 2) and inverted the relations in order 
to get a relation between the stellar mass to light ratio 
and the R band luminosity: 



M 

log(— ) = 7 log(LR) -f <^ 
Lr 



(7) 



for z = 0, we find: 



7 = 0.316 ±0.03 & <r = -2.7 



and for z = 0.8: 



= -13.1 



7 = 1.18 ±0.08 

We plot on Fig. [3] the lines corresponding to the range 
allowed by these relations. Obviously the total mass to 



light ratios Mtot/L derived in this work are systemati- 
cally higher than stellar M/L's inferred observationally. 
This can be understood by the fact that we consider 
the total mass, corresponding to both the baryonic and 
the dark matter components, whereas the observational 
AI/L ratio is the stellar mass to light ratio. 

We did the same exercise for the whole galaxy popu- 
lation (i.e. for all cluster-centric radii) to compare the 
{ri/2 , L) and (Mtot, L) trends inferred for the whole 
population to the one inferred for the central population 
(R< 500 kpc). When considering the whole galaxy pop- 
ulation, we find increasing (ri/2 , L) and (Mtot, L) trends 
at z = 0, both for COMA and viRGO . On the other hand, 
at z = 0.7, these trends are hardly defined and the scatter 
dominates, specially in the case of the (Mtot, L) trend. 

This comparison suggests that in the core of the cluster 
at high redshift, the galaxy population is already well re- 
laxed and constitutes a homogeneous galaxy population, 
whereas the galaxies at larger distances from the cluster 
center have more complicated dynamics. It appears that 
clusters have a 'well behaved relaxed nucleus' at high 
redshift. With time this relaxed behavior propagates it- 
self outward into the outskirts of the cluster, resulting at 
2; = in a overall 'well behaved' galaxy population. 



5. DISTRIBUTION OF STARS AND DARK MATTER IN 

INDIVIDUAL GALAXIES 

In order to get insight on how tidal stripping modifies 
the different galaxy mass components, we study the ratio 
of the total mass to the stellar mass and its evolution with 
redshift and cluster-centric distance. First, we verify that 
there is no dependence of the stellar mass with R . We 
find that the scatter in the stellar mass of the galaxy 
sample is large, and there is no discernable trend with 
Rat any redshift, suggesting that the stellar component, 
on average, remains roughly unaffected by tides. 

Then, we study the ratio of the total to the stel- 
lar mass, Mtot/M* and its evolution with both cluster- 
centric distance and redshift. Fig. 2] shows the result 
for z = and z = 0.7. The ratio of the total to the stel- 
lar mass clearly increases with cluster-centric distance. 
Sin ce tidal forces remove subhalo mass from the outside 
in (jPiemand et al.l [2007D . this result suggests that the 
extent of the dark matter subhalo (quantified by ri/2 ) is 
larger than the extent of the stellar component (quanti- 
fied by r*^2 , the half light radius) even in the very central 
part of the cluster. To investigate this further, we com- 
pare the extent of the stellar component r*^^ with rx/2 • 
We find that does not depend on the cluster-centric 
distance Rfor any redshift, both for COMA and VIRGO. 
In addition, for viRGO , we find that r^^^ is compact < 
2 kpc, whereas ri/2 > 8 kpc (Fig.[2l upper left panel, cen- 
tral bin). For COMA, we find that < 3 kpc, whereas 

ri/2> 11 kpc (Fig. O upper right panel, central bin). 
This shows that for even the most stripped halos in the 
central bin, ri/2 is still 4 times larger than r*^^ , sug- 
gesting that the dissipational stellar component is more 
compact and not as affected by tides. 

6. COMPARISON WITH GALAXY-GALAXY LENSING 

RESULTS 
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TABLE 1 

Deriving the scaling relations defined in Eq. [6] 



Redshift 









0.2 




0.4 




0.7 




MEAN 


COMA 


a = 
5 = 


0.315 ±0.111 
1.431 ±0.119 


a 
5 


= 0.165 ±0.136 
= 1.297 ±0.148 


a 
S 


= 0.274 ±0.090 
= 1.189 ±0.098 


a 
5 


= 0.194 ±0.059 
= 0.983 ± 0.077 


a - 


= 0.237 ± 0.045 
= 1.22 ±0.055 


VIRGO 


a = 
S = 


0.404 ± 0.050 
1.303 ±0.072 


a 
5 


= 0.213 ±0.036 
= 0.905 ± 0.068 


a 
5 


= 0.275 ±0.036 
= 0.963 ±0.064 


a 
5 


= 0.382 ±0.031 
= 1.186 ±0.053 


a 
S 


= 0.318 ± 0.02 
= 1.09 ±0.03 



In this Section, we compare the results from our simu- 
lated galaxy population with the results derived from ob- 
servational galaxy-galaxy lensing studies in clusters. We 
first need to be sure we compare things that are compa- 
rable. The quantity we use from the simulation is the 
3D half mass radius ri /2 ■ Former galaxy-galaxy lensing 
analyses performed through clusters have used a para- 
metric model to fit the lensing data, the so-called Pseudo 
Isothermal Elliptical Mass Distribution (piemd, see e.g. 
iLimousin et al.ll2005D . This profile is parametrized using 
two characteristic radii, the core radius and the scale 
radius (respectively rcorc and rcut in former articles). 
In galaxy-galaxy lensing analyses, since we cannot con- 
strain it, we usually fix the core radius to an arbitrary 
small value, making t his PIEMD profile equivalent to the 
IBrainerd et al.l (|1996( ) profile. As a result of the galaxy- 
galaxy lensing fit, we get some cons traints o n this scale 
radius. It is easy to show ( Eli'asdo ttir et al.|[2007.) that 
the scale radius for the piemd profile without a core ra- 
dius is equal to the radius within which the 3D mass 
equals half of the total mass. Thus this parameter in- 
ferred from the galaxy-galaxy lensing analysis can be re- 
liably com pared to ri /2 from the simulation. It can also 
be shown (|EHasd6ttir et al.ll2007| ) that for a PIEMD pro- 
file without a core radius, the radius within which the 
2D mass equals half of the total mass is found to be at 
3/4 of the scale radius. 

Of course, lensing is sensitive to the 2D projected sur- 
face mass density, and it may be confusing to state that 
we infer a scale radius related to the 3D properties of the 
halo from a lensing analysis. However, the 2D projected 
mass density we are sensitive to with lensing is a function 
of a scale radius that turns out to be the radius within 
which the 3D mass equals half of the total mass. 

We perform the comparison only with COMA since it 
is a 6 keV cluster and therefore more similar to the 
massive, lensing clusters studied observationally than 
VIRGO . Note that if need be, it is possible to rescale 
the r values found for the COMA galaxies to the tem- 
perature of a cluster probed observationally using the 
scaling relation: ri/2 0c T~^/^. The temperature of a 
cluster characterized by a velocity dispersion a scales as: 
T (X a^, and its mass M at a cluster-centric distance Tc 
scales as: M oc rc. From Eq. [2] and [21 we can derive 
that the tidal radius of a galaxy orbiting in the cluster 
with a circular velocity Vc scales as: 

2/3 1/3 

so Vf OC VcTcl cr, and finally oc T^^l"^. This means 
that the more massive a galaxy cluster is (thus the more 
dense the environment), the more severe the tidal trun- 



cation. For the purpose of the comparison with observed 
clusters, all apparent magnitudes have been converted to 
the absolute rest frame R band magnitude. We will ap- 
ply the same selection criteria (magnitude cut-offs) to the 
simulated galaxies in order to compare the same luminos- 
ity galaxy population to the one probed in galaxy-galaxy 
lensing studies. Finally, the galaxy-galaxy lensing results 
have been rescaled to the median value of the luminos- 
ity of the simulated galaxy population, using the scaling 
laws defined in Eq. [H] with a = 0.5 and 8 = 1.0. 

In particular, we compare the resul ts presented in 
the pre vious sections with the work by Limousin et al.l 
()2007al) . which probes the cluster population down to 
a magnitude of -17.5 for a sample of five galaxy clus- 
ters at redshift around 0.2, with a mean X-ray tempera- 
ture of 8.5 keV. As discussed in detail in ILimousin et ahl 
(|2007al) . these detections are in agreement with the sim- 
ilar studies performed by Natarajan et al. on observed 
clusters r anging from z = 0.17 - 0.58 with space based 
HST data (|Nataraian et al.|[2007[) . Fig.[5]shows the com- 
parison. We find qualitative agreement. However, due 
to the weakness of the galaxy-galaxy lensing signal and 
the limitations of ground based data, error bars on the 
estimated parameters are still rather large. Also shown 
in Fig. [5] is the slope derived for the trend of the vari- 
ation of typical subhalo mass (mass of a subhalo that 
hosts an L* early-type galaxy) as a function of cluster- 
centric radius for the massive lensing cluster CI 0024+16 
derived from space-based data. The data-set comprises 
of an HST mosaic that samples out to 5 Mpc from the 
cluster center and the constraints on the subhalo mass 
are derived applying galaxy-galaxy lensing methods in 
3 radial bins. Strong and weak lensing constraints are 
combined to derive the average properties of a typical 
dark matter subhalo that hosts an L* early-type galaxy 
in each radial bin. The variation of the mass of the fidu- 
cial subhalo with cluster-centric radius is well-fit by a 
linear relation. While the slope of this relation can be 
directly compared to that derived from the simulated 
COMA cluster (as shown by the dashed line in Fig. 4) the 
normalization cannot be compared as the central den- 
sity of CI 0024+16 is significantly higher. We find that 
the trend derived from these space-based lensing observa- 
tions is in agreement with the results of the simulations. 

7. CONCLUSIONS 

Using high resolution, hydrodynamical N-body simu- 
lations of two fiducial galaxy clusters, one with param- 
eters typical of COMA and the other viRGO, we study 
the tidal stripping process in detail. These are the first 
simulations that include the baryonic component that 
have been used to study this process, previous studies 
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Fig. 3. — Upper panels: median half mass radius r^j^ as a function of luminosity for the central galaxies sample (R < 500 kpc), for 2 = 
(solid squares )and 2 = 0.7 (dashed circles); VIRGO (left) and COMA (right). Shown is the best fitting linear relation to the unbinned data 
points (Table [ij. We find that ii^creases with luminosity, but has a larger scatter in the higher redshift slice. Lower panels: median 
total mass to light ratio as a function of luminosity for the central galaxies sample f or 2 = and z = 0.7; V IRGO (left) and COMA (right). 
Overplotted lines correspond to the range allowed by the fundamental plane analysis bv lJergensen et al.l l|2006l) . Note different axis scalings 
in the figure panels. 
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were limited to dark matter only simulations. We in- 
fer a strong trend between the extent of cluster galaxy 
dark matter subhalos quantified by the half mass ra- 
dius ri/2and the projected cluster-centric distance R. 
We show that the dark matter component is preferen- 
tially stripped, whereas the stellar component is much 
less affected by tidal forces. We infer a trend in these 
simulations that is much stronger than the one inferred 
from the d ark matter o nly co smological N-body simula- 
tions of Ghig na et al.l ()1998D . This could suggest that 
tidal stripping is more efficient in the inner regions of 
clusters when the effects of baryons are included. How- 
ever, we caution that comparing both simulations is not 
easy and therefore we cannot assess reliably the impact 
of the baryonic component in the tidal stripping process 
at this point. What is needed now is a devoted dark 
matter only simulation of viRGO and COMA in order to 
compare reliably the outputs of each study and investi- 
gate the expected influence of the baryonic component 
on the tidal stripping process. 

With the next generation of space telescopes, in partic- 
ular with wide field space based imagers, such as dune"* 
and SNAP^, we will be able to probe this trend in a 
sample of galaxy clusters spanning different range of 
masses and dynamical states. The limiting aspect of 
these kind of studies are the relatively small mass of 
the dark matter subhalos associated with cluster galax- 
ies (leading to a small modulation in the overall tan- 
gential shear field making detection challenging), and 
therefore the large number of background sources with 

* www.dune-mission.net 
5 www.snap.lbl.gov 



reliable shape measurements needed for the lensing anal- 
ysis. A very promising technique will be the use of fu- 
ture radio and millimeter-wave interferometers: accurate 
measurements of the detailed dynamical structure of the 
background galaxies, in particular rotating disks, should 
make it possible to probe the shear directly on each in- 
dividual galaxy, thus allowing shear measur ements with 
much less gala xies than it is currently done (|Blainll2002t 
IMoralesll2006l) . 

From the numerical point of view, what is needed is 
a sample of galaxy clusters that span a wide range in 
mass and dynamical states in a large simulation box in- 
cluding baryonic physics, so that a robust comparison 
between observations and theoretical expectations can 
then be made. Such a quantitative comparison will be 
possible in the near future and will provide further in- 
sights into the physics of cluster assembly and the process 
of tidal stripping. 
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Fig. 4. — Median ratio of the total mass to the stellar mass of galaxies as a function of projected cluster-centric distance Rat 2; = and 
z = 0.7. Left: COMA ; right: VIRGO . 
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Fig. 5. — Comparison of Limousin et al. (2007a) and Natarajan et al.(2007) galaxy-galaxy lensing results with simulated COMA galaxies. 
Left panel: versus cluster-centric radius; Right panel: total mass versus cluster-centric radius. To enable direct comparison with the 
lensing data, the same magnitude cut-off has been adopted in the simulations. The data point with error bars shown is derived from galaxy- 
galaxy lensing results averaged over the subhalo population ranging in cluster-centric radius from to 1.8 JVIpc for a sample of massive 
lensing clusters. The slope indicated with a dashed line shows the trend derived for CI 0024-1-16 from galaxy-galaxy lensing techniques as 
well by Natarajan et al. (2007) using space-based data. This trend is derived from mosaiccd HST WFPC-2 data for CI 0024+16 out to 5 
Mpc from the center. The slope corresponds to the change in subhalo mass with cluster-centric radius for a typical subhalo that hosts an 
early-type L* galaxy. Note that the normalization cannot be compared directly due to the fact that the central density of CI 0024-1-16 is 
significantly larger than that of the simulated COMA studied here. 
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